Cosmic Axion* 



Jihn E. Kim 

Lyman Laboratory of Physics, Harvard University, Cambridge, MA 02138, and 
Department of Physics, Seoul National University, Seoul 151-742, Korea ^ 

Abstract 

I review the axionic solution of the strong CP problem and current status of 
the cosmic axion search. 

I. INTRODUCTION 

Quantum chromo dynamics before 1975 considered the following Lagrangian 

C = --^TrF^.F^'^ + q{tD^7^ - M)q. (1) 

where M is the diagonal, 75-free, real quark mass matrix. But after 1975, the following term 
is known to be prensent in general in a world without a massless quark, 

+ Y^TrF,.F- (2) 

Since this 9 term violates CP invariance, the upper bound of the neutron electric dipole 
moment puts a strong constraint on the magnitude of 9, \9\ < 10~^. The smallness of 9 has 
led to the strong CP problem, "Why is 9 so small?" 0] 

We know that many small parameters in physics have led to new ideas, in most cases 
leading to new symmetries. For example. My/ /Mp <^ 1 has led to supersymmetry, niu^ <C 1 
GeV has led to SU{2)l x SU{2)r chiral symmetry, etc. For the strong CP problem, the 
nicest solution is the very light axion resulting from the Peccei-Quinn symmetry p|. 
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II. THE AXION SOLUTION 



The reason that the axion solves the strong CP problem is the following. This argument 
is due to Ref. 0. In the axion solution, ^ is a dynamical field, but for a moment let us treat 
it as a parameter (or coupling constant). The partition function in the Euclidian space after 
integrating out the quark fields is 

^-y\S] ^ J idA^] Det(D'^7M + m.) ^M' J ^H^^F' - %Q{FF})\ (3) 

where { } includes the factor l/327r^. It is known that Det factor in the above equation 
is positive Also note that the d term is pure imaginary. Therefore, using Schwarz 
inequahty, we obtain the following inequality 

< Det(Z}'^7M + ^i) e^ph / d^xi^F-" - %Q{FF})\ 

= J[dA^] U^ Det(D'^7^ + m,) exp[- / d^x^F'] = g-Z'^'^^M (4) 

which gives 

V[9] > V[0]. (5) 

Thus ^ = is the minimum. However, if ^ is a coupling constant, any 6 can be a good 
coupling constant as any «em is allowed theoretically. The axion solution interprets ^ as a 
dynamical field, introducing a kinetic energy term for for the boson field 6. In this case, we 
necessarily introduce a mass parameter Fa, accompanying the axion a, 

e-y. (6) 

Then the shape of the potential of a is as shown in Fig. 1. 

The hight of the potential is guessed to be of order Aq^j-,. The current algebra calculation 
gives {2Z/ {\ + Zy')flm^. Since the instanton solution gives / d!^x{FF} = Z and FF appears 
in the form given in Eq. (2), 6 is a periodic variable with period 2tt. Since ^ is a dynamical 
field, different ^'s do not describe different theories, but merely different vacua. Thus, as 
universe evolves, 6 seeks the minimum of the potential ^ = 0. This mechanism explains 
very elegantly why 6 is so small in our universe. The above proof assumed no CP violation 
except that from the 6 term, and the weak CP violation shifts the minimum point very 
little, 6 ~ 10^^'' [0 which is far below than the bound given by the neutron electric dipole 
moment. 

An important feature is that a does not have any potential except that coming from 
6{FF}, otherwise the mechanism does not work. The effect of weak CP violation introduces 
a potential, but as commented above the effect is very small. 

To make 6 dynamical, one must have a mechanism to introduce a scale Fa. Depending 
on the nature of a, one can classify axion models into three broad categories: 

(i) a is the Goldstone boson of a spontaneously broken chiral U{1) symmetry. The divergence 
of this U{1) current must carry the color anomaly oc FF so that {a/Fa){FF} coupling 
arises. 
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FIG. 1. A schematic view of the axionic potential with a minimum at 6 = 0. 

(ii) a is a fundamental field in string models. The scale Fa arises from the compactification. 
It is called the model-independent axion 

(iii) a is a composite field. Fa arises at the confinement scale [^. 



A. Domain walls 

Because ^ is a periodic variable, the axion potential looks like as the one shown in Fig. 2. 
In this example, the origin < a >= is identified with the vacuum < a >= GvrFa both of 
which are denoted as black dots. Thus, < a >= 0, 27: Fa, and An Fa are the three degenerate 
vacua, distinguished by a black dot, a star, and a triangle. Since the discrete symmetry of 
vacua is spontaneously broken in the evolving universe, there appear three kinds of domain 
walls, i.e. Njj^r = 3 in our example, in the evolving universe. This leads to the so-called 
axionic domain wall problem However, if Ndw = 1, there seems to be no domain wall 
problem even if they are formed in the evolving universe. This is because the string domain 
wall network system in the Ndw = 1 model can be erased easily. A large string attached 
with a large domain wall dies out due to punched holes in the wall expands with light velocity 
erasing the wall. 

If a singlet scalar field develops a VEV v, usually the axion coupling to the gluons has 
the form 

(v/Nnw) ^ ^ 

which implies that the coupling is smaller by a factor N^w- Thus the axion mass is larger 
by a factor N^w if one uses the vacuum expectation value of the Higgs field. However, if 
one uses F^, there does not appear the dependence on Ndw as is evident from the definition 
of Fa in Eq. (6). One can imagine a possibility that 10^^~^^ GeV scalar vacuum expectation 
value with Ndw ~ 100 can be consistent with the cosmological bound. But in this case, of 
course. Fa ~ lO^^ q^v. 
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FIG. 2. An axionic potential with Now = 3. Here 9 = a/ Fa- 



B. Superstring axion 

String models include massless bosons Gmn {MN = symmetric), Smat {MN = 
antisymmetric), and dilaton 0. Among these, Bmn is of our interest here. Any D di- 
mensional index M can take {D — 2) transverse directions for a massless particle. Therefore, 
Bmn has {D — 2){D — 3)/2 physical degrees. In 4 dimensional Minkowski space time, 
B^i, {fi, u = 0,1, 2, 3} has one physical degree; thus it is a pseudoscalar. The pseudoscalar 
is the dual of the field strength of B^^, 

d.aui ~ t^upaH^"' (8) 

where aMi is called the model- independent axion (Mia) in string models @. The Mia 
coupling is universal to all fermions 

H^,pi,^^'^'Yi^ ~ d^aMii'Yl^i'. (9) 

Of course, the coupling of qmi is only of the derivative form, rendering the nonlinear global 
symmetry, a a+ (costant). This symmetry is anomalous and the Mia coupling is universal 
to all gauge groups [| 

~ a{FF + F'F' + ■■■) (10) 

Since any superstring model possesses Mia, the axion solution of the strong CP problem 
gets a firm theoretical support in string models. But the axion decay constant is too big 
in a naive string models [§. In anomalous U{1) models, however, the axion decay constant 
can be lowered 



^This comes from the relation H = dB +uj3y — ^3L where losy and lo^l are Yang-Mills and Lorentz 
Chern-Simons three forms. 
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III. AXION PROPERTIES 



Remembering that the axion is a dynamical 9, we can easily derive its interaction terms. 
For this, we follow a simple route of effective field theory. 

The simplest axion example is the heavy quark axion |TD|. Note that the axion models 



should provide a pseudoscalar a, coupling to FF. The a is housed in the complex scalar 
singlet field a. By introducing a heavy quark Q, the following Yukawa coupling is introduced, 

C(xaQRQL + \i.c. (11) 

This model posseses a global Peccei-Quinn symmetry, Ql — > c'^^/^Ql, Qr e'^/^Q/?, o" — >■ 
e*°(T, and 9 9 — a. The VEV < a >= Fa/ gives a mass to Q, and produces a Goldstone 
boson a where a = {l/^/2){Fa + p)e^°'^^°'. Below the scale Fa, the light fields are the gluons 
and a. The Lagrangian respecting the above symmetry is 

C = Ud^af + (derivative terms of a) + ^{9 + ^)F;^F''^r (12) 

Thus, minimally we created a dynamical variable 9 = 9 + {a/ Fa). It is redefined as a/ Fa by 
shifting the a field. From now on, 9 implies 9. 

Next, let us introduce the known light quarks. As the first extension, let us consider the 
up quark condensation in one-flavor QCD. The mass term in this theory is given by 

Craass = -niuURUL + h.C. (13) 

Formally, we can assign the following U{1) chiral transformation. 



(14) 



ul — 








Ur — 


- e-^'uR 






m — 




9 - 


^9 + 2a 



Due to the above chiral symmetry, we expect the following effective potential below the 
chiral symmetry breaking scale ^qcd, 

+ A4-r#-e2^^-2^^ + --- + h.c. (15) 

where ■ ■ ■ is the higher order terms, A's are couplings of order 1, (uu) = f^e*'"'^, and the 
QCD scale lyqcD is inserted to make up the correct dimension. In addition, 6=*=*^, e^*"^*^, etc 
is multiplied to respect the U{1) symmetry. Note that if m„ 7^ and 9 is not a dynamical 
variable, then the strong CP problem is not solved. Note that, if = then only the 
rriu-independent terms survive, leading to 

V = -lAiAt;3e^S-^ + \^t-e^^'--^'' + ■■■ + h.c. (16) 

Thus, redefining the r] field as rj' 
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f]' = 1] — v6, 



(17) 



the 9 dependence is completely removed from V. The 6 parameter is unphysical if a quark is 
massless. Namely, the massless up quark scenario solves the strong CP problem even though 
it obtains a constituent quark mass. The relevance of this solution hinges on the viability 
in hadron physics phenomenology [|TT|. For m„ 7^ 0, at the minimum < a >=< rj >= 0, the 
mass matrix is 



( XKqcdv + X'muv -hh3£^ \ 

V Fa F2 



(18) 



It is easy to calculate determinant of M 
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DetM^ = !%i^(AA').^ - XAlcD - >^'mK%cD- (19) 



For Fa ^ others, we obtain = [XAqcd + X'm)v. Thus the axion mass is 

- A^cdI (20) 



2 _ rriuAQCD ( XX'v"^ 2 

F2 {xAqcdv + X'tu^v 



which is supposed to be positive. Otherwise we should have chosen a = irFa- This axion 
mass shows the essential feature: it is suppressed by Fa and multiplied by m„. The rest 
is the condensation parameters. Usually, the condensation parameters are given in two or 
three quark flavors. 



A. The invisible axion mass 

For two flavors of u and d, we can repeat the above argument with U{l)u x f/(l)d 
symmetry 

ul e"^UL, ur e~"^UR, di e'^dL, d^ e'^-f^dn 

ruu e-^'"mu, ma e'^^^mrf, 9 ^ 9 + 2{a + /3). (21) 

The effective potential respecting the above symmetry is 

V = \mumdA^cD^^^ ~ — ^ "^"^ ^ ^ ~ \X'2fnu < uu > — ^Agm^ < dd > 

- |A^m; <dd> e-'^ - iA'gm* < uu > e"^^ + • • • + h.c. (22) 

We can diagonalize the 3x3 mass matrix of a, tt^ and rj where the phases of < uu > and 
< dd > are proportional to ^ + ^ and ■§ — respectively. 9 is proportional to The 
axion mass is 
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where Z = rriu/md. The above mass formula is vahd for the very hght (or invisible) axion. 
For the PQWW axion we need an extra consideration of separating out the longitudinal 
degree of the Z boson. Below the chiral symmetry breaking scale the axion Lagrangian is 

£ = -(<9^a)^ — TJ^a^ + (interaction terms). (24) 
The interaction terms depend on models. 



B. The KSVZ axion 



The KSVZ axion [10| introduces a heavy quark Q, 



- £ = faQuQi + h.c. + rriuURUL + madadL H (25) 

where a provides a. The light quarks does not transform under the shift of a. At tree level, 
there does not exists an axion-electron coupling and it can be induced at one-loop order. 



C. The DFSZ axion 



The DFSZ axion [0 introduces two Higgs doublets 

- £ = XaaHlH; + fuURULH^ + fdndLHl + /^e^e^i/O + h.c. (26) 

where a resides mostly in a with a small leakage to and ; phases of {H^-, ~ 
{cos P, sin P} a/ Fa where tan/3 = ^2/^1- Depending on models, Hi,H2, or the third Higgs 
doublet Hs can couple to the electron. For the first two cases, the electron coupling arises 
at tree level, ~ {cos /3, or sin P} {a/ Fa)meei'j5e. 



D. The a — 7 — 7 coupling 



In view of the possible detection of the cosmic axions in a high-q cavities immersed in 
the strong magnetic fields [|13|, it is important to know the axion-photon-photon couplings. 



More than a decade ago [|1J], it was calculated, but the current citation of the coupling is 
not accurate. The details of the KSVZ and DFSZ couplings are given in Ref. [l^. The 
chiral symmetry breaking at 100 MeV shifts the 077 coupling. Thus the coupling is usually 
expressed as 



2 A + Z 

3 ' 1 + Z 



c, 



1.92 



(27) 



where we used Z 



mu/rud 
. The c 

given by the Peccei-Quinn symmetry of the theory, 



The Ca-y-y is the coefficient of 



{a/Fa){FemFem} term. The c, 



0.6 in the last equation 
a-y-y is given above the chiral symmetry breaking scale, and is 
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E 



Ca77 = where E = TiQ^^Qpq, 5abC = TrXaXbQpQ- 



(28) 



The normahzation is such that the index / for 3 and 3* is |. The Peccei-Quinn charge is 
derived from the currents obtained from the Lagrangians given in Eqs. (25) and (26), 



KSVZ : J^Q = d - iQ7M75Q 

DFSZ : JPQ = v + J2i Ui-if.lbUi + Ei dat^^bdi 



(29) 
(30) 



where a; = ^2/^1 is the ratio of the Higgs doublet VEV's. It is given in Table 1 [IT3. Here 
cr denotes the electric charge of the representation R in units of the positron charge. 

Table 1. The axion-photon-photon couplings. 



KSVZ 




DFSZ 




eR 


Ca-yy 




Ca-yy 


eR = 


-1.92 


any {i = 1) 


0.75 


63 = -| 


-1.25 


1 (z = 2) 


-2.17 


63 = 1 


0.75 


1.5 {i = 2) 


-2.56 


63 = 1 


4.08 


60 {i = 2) 


-3.17 


68 = 1 


0.75 


1 (z = 3) 


-0.25 


"^3 3)3 


-0.25 


1.5 {i = 3) 


-0.64 






60 {i = 3) 


-1.25 



The above table cites the couplings in the KSVZ and DFSZ toy models. In reality, there 
can be many heavy quarks which carry nontrivial Peccei-Quinn charges, e.g. as in Ref. 0. 
For example, superstring models usually have more than 400 chiral fields. Also, the light 
quarks are most likely to carry the Peccei-Quinn charges. Therefore, these effects add up. 
In superstring, different models give different values for Ca-y-y- If the standard string model 
is known, we can predict the exact value of c^^^ in such a model. 



IV. ASTROPHYSICAL BOUNDS 

For a sufficiently large Fa, axions produced in the stellar core escape the star easily, 
which provides an efficient way of the energy loss mechanism. Comparing this axion emission 
process with the standard energy loss mechanism through neutrino emission gives a bound 
on Fa- The production cross section is the dominant bottle neck. Thus the stellar bound is 
such that enough axions are not produced, giving a lower bound on Fa- Any axion model 
has the Primakoff process of the axion production and nucleon collision process of the axion 
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production. The DFSZ model has additional Compton type axion production and similar 
electron (or positron) scattering axion production processes. But for Fa > 10^ GeV which 
is of our interest here, both the KSVZ and DFSZ models have similar lower bounds. The 
astrophysical bounds are reviewed extensively in the literature |jl6|. 
one obtains F„ > 2.3 x 10^ 



For example, from Sun 
(ca77/0/75) GeV, from red giants one obtains 0.9 x 10^ GeV, 



from globular clusters one obtains 2 x lO^Ca^^ GeV. For the supernova, Iwamoto and others 



studied before the discovery of SN1987A [|T7|. But these pre-SN1987A papers failed to give 
a strong lower bound. After SN1987A, many groups obtained the lower bound of order 
10^- 



-10 



GeV 18 



The discrepancy of the numerical studies before and after the discovery of 
SN1987A was due to the axion couplings used in the analyses. Of course, the correct coupling 
is of the derivative form duaNj^'j^N with nucleon A^. The correct Goldstone boson nature 

In general, this consideration of 
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of the pion is also important as pointed out in Ref. 
the derivative coupling is important at high temperature as in the supernovae, and gives a 
lower bound of order 



Fa > 10^ GeV. 



(31) 



V. COSMIC AXION 

The U{1) global symmetry breaking is achieved by a Higgs potential shown in Fig. 3. 
The circle at | < cr > | = is the axion oscillation direction. The small perturbation at 
\ < a > \ = Fa arises due to QCD instanton effects. In Fig. 3, there are three degenerate 
minima, leading to an N^w = 3 model. In the evolving universe, the U{1) breaking at Fa 
produces axionic strings. At the chiral symmetry breaking scale of ~ 100 MeV, the domain 
walls are attached to these axionic strings, which is shown in Fig. 4. The axion-string and 
axion-domain-wall system does not die out quickly if N^iy ^ 1. 




FIG. 3. A potential breaking C/(1)pq with Ny)^/^/ = 3 arising at the QCD phase transition. 
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FIG. 4. A schematic view of cosmic string and domain walls with Njjiy = 3. 



For some time, the reheating temperature after inflation Trh is greater than Fa so that 
the baryon number generation through GUT interactions dominates. In this case, we must 
allow only models with A^£)vi/=1- 

But the condition for Now = 1 is not necessary if the reheating temperature after 
inflation falls below Fa- In supergravity models, the gravitinos which interacts very weakly 
with observable sector particles decay so late in cosmic time scale that they can dissociate the 
preciously nucleosynthesized light elements unless Trh < 10^~^° GeV [^. These arguments 
favor a low reheating temperature, presumably below Fa- Then the axionic strings and 
domain walls are not important at present. 

Thus the reliable constraint from cosmology is the cold axion energy density ||22|. This 



arises from the reason of the invisible axion's extremely feeble coupling so that its lifetime 
is many orders larger than the age of the universe. The almost flat axion potential is felt 
in the evolving universe when the Hubble parameter becomes smaller than the axion mass, 
3H < rUa- This condition is satisfied at the cosmic temperature ~ 1 GeV. In the inflationary 
universe, the vacuum value of < a > is the same in the visible universe. At T ~ 1 GeV, 
< a > begins to roll down the potential hill and continues the oscillation around < a >= 0. 
This motion of the coherent axion oscillation carries a huge energy density and behaves 
like nonrelativistic particles. By now, the estimate of these cold axion energy density is 
standard, and one obtains [^1 




1.18 



Qah' ^ 0.13 X W^'-'A^^'m) ^^^^ Nl^ (32) 



where 



,„„ = » " ' f eV. (33) 



F 
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A200 is the strong interaction scale in units of 200 MeV, Qa is the axionic fraction in critical 
energy density, h is in units of 100 km/s/Mpc, 61 is the initial misalignment angle taken 
as 01 ~ vt/a/S, and f{9i) is a correction factor of order 1. The above consideration gives 
F'a ^ 10^^ GeV not to close the universe by the cold axions. 

If = 6, then the axion mass to close the universe is 0.8 x 10~^ eV. 

Superstring models give Now = 1- For Qa ~ 0.7, h ~ 0.65, the axion mass is 0.5 x 10~^ 

eV. 

Let us mention the axion string and domain walls in standard Big Bang or inflation 
with Trh > Fa- In this case, for N^w > 1 a complicated axionic string and domain wall 
network do not die out easily. Therefore, only Ndw = 1 models are viable. Even for the 
Now = 1 case, the string-wall system generates considerable energy. One group asserts 
that the string- wall system outweighs the cold axions p4|, while the other group advocates 



energy density of axions from walls is comparable to or smaller than the cold misalignment 



axions 



The difference comes from the different assumptions on the nature of axions 

Recent estimate of the 



radiated from the axionic walls whether they are cold p4| or hot 
axion energy density from axion walls by Battye and Shellard is dominated by the axionic 
string loops [El], 



^string loop ^ IQ.T/l^^A 



'1 + ")3/2 
K 



Fa 



1012 GeV 



1.18 



(34) 



where a is roughly denoting the loop creation size relative to the horizon size and k is the 
back reaction decay rate. For a ~ k ~ 0.1, we have rria > 100 fieV or < 4 x 10^° GeV. 
On the other hand, a few years ago Nagasawa and Kawasaki gave a stronger bound 
Fa < 10^° GeV which results from large strings domination over the loops. In any case, with 
Trh < Fa in inflationary models, this string and wall consideration is irrelevant. 



VI. AXION SEARCHES 

The axion search is really on the invisible axion closing the universe, for Fa ~ 10^^ GeV. 
In this case, the axion mass to be searched for falls in the several fieV region. For models 
with = 6, the vacuum expectation value of the singlet Higgs field is around lO^^^ GeV 
region. 

The basic assumption for the cosmic axion detection is that the cold axions comprise 
about 70% of the closure mass density of the universe. In our galaxy, it is about 0.3 GeV/ cm^. 
Even though the axion interaction is very weak, the enormous number density overcome the 
extremely small conversion rate of axions to photons. This cosmic axion detection employs 
Sikivie's high-q cavities. (Note, however, the Univ. of Tokyo effort to search for nuclear Ml 



transitions ||2^ which does not employ Sikivie's detector.) The photons produced by the 
cosmic axion Primakoff process in the strong magnetic field are collected in the cavity. There 
already exist the first round experiments (the Rochester-Brookhaven-Fermilab (RBF) group 
and the University of Florida (UF) group experiments tried to detect photons collected in 



this cavity [^]. They are shown at the upper right-hand corner in Fig. 5. The current 



experiment at LLNL repeats the same type of experiment but with a bigger cavity [29|. The 



sensitivity of this new LLNL experiment as of June, 1997 is also shown in Fig. 5. Next year, 
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the sensitivity of the LLNL group increases considerably as shown with a bigger box in Fig. 
5. 

On the other hand, the Kyoto group passes the Rydberg atoms in another cavity where 
the photons collected in Sikivie's cavity are shone into the Rydberg atoms; the electrons 
freed from the Rydberg atoms by the axion-converted photons are measured PU . 



In Fig. 5, we also show a few model predictions in the KSVZ and DFSZ models. The 
DFSZ in Fig. 5 represents the {df^, e) unification model. 

In Fig. 6, we compare these data with the predictions from several very light axion 
models. As is evident from the figure, it will be difficult to pin point a toy model even 
if the cosmic axion is detected. Most probably, the very light axion would come from the 
Pecce-Quinn symmetry breaking where both heavy quarks are the light quarks carry Peccei- 
Quinn charges. Superstring models have this property 0. If a standard superstring model 
is known in the future, the axion detection rate will be predicted in the axion dominated 
universe without ambiguity. 

The detection of the cosmic (or galactic) axions would be a stunning confirmation of both 
the particle physics theory (instantons, invisible axions, etc.) and modern cosmology (galaxy 
formation, dark matter, etc.), and would open an window toward a possible fundamental 
theory of everything. 
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FIG. 6. A part of Fig. 5 compared with several invisible axion models. 
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